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Large amplitude collective motion is investigated for a model pairing Hamilto- 
VO ' nian containing an avoided level crossing. A classical theory of collective mo- 

tion for the adiabatic limit is applied utilising either a time-dependent mean- 
field theory or a direct parametrisation of the time-dependent Schrodinger 
, equation. A modified local harmonic equation is formulated to take account 

\ of the Nambu-Goldstone mode. It turns out that in some cases the system 

selects a diabatic path. Requantizing the collective Hamiltonian, a reasonable 
agreement with an exact calculation for the low-lying levels are obtained for 



C***"" ■ both weak and strong pairing force. This improves on results of the conven- 

tional Born-Oppenheimer approximation. 
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I. INTRODUCTION 



Nuclei are finite fermionic many-body systems which support many kinds of collective 
motion. While properties of high-frequency vibrations (giant resonances) may be well re- 
produced in the small-amplitude limit by the random-phase approximation (RPA), some 
low-frequency vibrations exhibit a strongly anharmonic nature that the RPA cannot de- 
scribe. Nuclear fission and shape coexistence phenomena also have such a large amplitude 
nature. In order to investigate these kinds of large amplitude collective motion, one would 
like to reduce the number of degrees of freedom to a few judiciously chosen slow collective 
coordinates. However, in nuclear systems, this task is not trivial, since there is no obvious 
separation of scales. 

In molecular physics, where the masses of electrons are so much smaller than those 
of the atomic nuclei, the electronic motion is normally much faster than nuclear motion. 
Thus, the collective coordinates are usually functions of nuclear coordinates and the Born- 
Oppenheimer (BO) approximation works very well. On the other hand, in nuclear physics, 
since a nucleus consists of neutrons and protons which have almost the same masses, both 
the definition of collective coordinates and the applicability of adiabatic assumption are 
never obvious. 

Although a large number of studies have been done to calculate the potential energy 
surface using the constrained Hartree-Fock (HF) or Hartree-Fock-Bogoliubov (HFB) theory 
with a given generalised cranking (or constraint) operator, the choice of collective coordinate 
(i.e., the choice of cranking operators) has been rather arbitrary in most cases. 

In this paper we shall apply a special theory Q, that is designed to determine a self- 
consistent cranking operator for adiabatic large amplitude collective motion (ALACM). This 
theory provides a method to find approximate decoupled motion which is confined to a 
few dimensional submanifold of the configuration space, within the framework of classical 
Hamiltonian dynamics. Since most systems of practical interest are not exactly separable, it 
is important that the theory can provide a submanifold which is approximately decoupled. 
This method has been applied to problems in nuclear physics @,[| and in other fields (see 
references in [[!]]) 

In nuclear collective motion which involves level crossings, the problem of adiabaticity 
versus diabaticity has been a long standing question since this problem was discussed by 
Hill and Wheeler |4]]. During a nuclear shape change, the diabatic process is often more 
favoured than the adiabatic one ||. This raises the question whether an adiabatic theory, 
such as ALACM, can be used to shed some light on diabatic dynamics. An answer to this 
question is one of main goals of this paper. 

In nuclear phenomena, it is well-known that the pairing (superfluidity) influences all low- 
frequency collective motion. A well-known example is the effect on the moment of inertia 
for rotational nuclei, which is always smaller than the rigid-body value at low spin, which 
can be explained as an effect of pairing correlations. At the same time, the ground states of 
heavy nuclei with open-shell configurations are reasonably well described by the superfluid 
Bardeen-Cooper-Schrieffer BCS wave functions with energy gaps of about 1 MeV. Properties 
of both collective and non-collective (quasiparticle) excitations depend on size of the energy 
gap. Furthermore, it has been argued by Bertsch || that nuclear shape change may be 
associated with the hopping of nucleon pairs by means of the pairing force. In this case, the 
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pairing interaction and level crossings play an essential role to determine the collective mass 
( "hopping mass" ) . 

According to these considerations, the pairing interaction should play a key role in un- 
derstanding the large amplitude collective motion in nuclei, especially when level crossings 
are involved as the shape change is taking place. Therefore, it is important to investigate 
the applicability of ALACM for such a system with level crossings and a pairing force. The 
theory has not been applied to such systems before. 

In this paper, we study a model Hamiltonian describing a system interacting through 
the pairing force. The model has a single-particle level crossing and multiple local minima, 
and thus may be regarded as a model for shape-coexistence phenomena. In Sec. |Tj, the 
formalism of ALACM is briefly recapitulated. In order to apply the ALACM to a classical 
Hamiltonian with a spurious component (Nambu-Goldstone mode), a modified version of 
the local harmonic equations is formulated. The classical Hamiltonian for a pairing Hamil- 
tonian is derived in Sec. [TTT| , both by using the mean-field (in this case BCS) theory and by 
applying a parametrisation which exactly conserves particle number. The results of numer- 
ical calculations for a simple two-level system are given in Sec. [IV] and the conclusions and 
an outlook are summarised in Sec. [Vj. 

II. BRIEF REVIEW OF ALACM 
A. Local harmonic equations (LHE) 

We briefly review the theory of ALACM (see Ref. Q for a complete description). In this 
section, we use a summation convention where the repeated appearance of the same symbols 
(a, (3, ■ ■ • ; i, j, • • •) in upper and lower indices denotes a sum over that symbol for all possible 
values. We also use the convention that a comma in a lower index indicates the derivative 
with respect to the coordinate, thus F >a = dF/d^ a . 

The theory of adiabatic large amplitude collective motion (ALACM) is applicable to a 
classical Hamiltonian system which has kinetic terms only quadratic in momentum. We 
thus have to start with a truncated Hamiltonian 

H(& tt) = + V(0 , a, (5 = 1, - ■ ■ , n , (2.1) 

where the mass tensor B alS , in general, depends on the coordinates £° and is defined by 
truncation of the Hamiltonian to second order 

dli 0,071/3 

Thus all terms more than quadratic in momentum are neglected. In the sense that the 
higher-order terms are small, this theory may be regarded as an adiabatic theory in the 
small- velocity limit. The tensor B a p, which is defined as the inverse of B af3 [B ai B 1 p = 8g), 
plays the role of metric tensor in the Riemannian formulation of local harmonic equations 
(LHE) below. 



(2.2) 

7T = 
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Collective coordinates <f and intrinsic (non-collective) coordinates q a which are approx- 
imately decoupled from each other, are assumed to be obtainable by making a point trans- 
formation, conserving the quadratic nature of Eq.(|2.1|). 



9* = /*(fl (i = l, •••,#). (2.3) 
q a = f a (0 (a = K+l,-.-,n). (2.4) 

In this section, we use symbols (a,/3, ■ ■ ■) for indices of original coordinates, (/x, u, ■ ■ •) for 
new coordinates after the transformation, ■ ■ ■) for collective coordinates and (a, b, ■ ■ -) 
for intrinsic coordinates. The new Hamiltonian after the point transformations takes the 
form, 

H « \&i VlVj + l -B ab PaPb + V{q\ q a ) , (2.5) 

where we have used the fact that for decoupled motion the mass tensor must be block- 
diagonal, 

B ai = . (2.6) 

Besides the block-diagonality of mass tensor, we require the absence of both "real" and 
"geometrical" forces orthogonal to the decoupled manifold, 

V a = , (2.7) 
B% = . (2.8) 



In practice, these three decoupling conditions, (|2.6|) , (|2.7|) and (|2.8| ), cannot be satisfied 
exactly, except for special cases. Thus, we need a method applicable to an approximately 
decoupled manifold. The Riemannian formulation of LHE ]I[ is the one we choose to use in 
this paper. In a case of a single collective coordinate (K =1), the basic equations of this 
formalism can be written as 

V* = A/i , (2.9) 

wi = ^/i ■ ( 2 - 10 ) 

Here the covariant derivative (denoted by ;) in the l.h.s. of Eq. ([2.10 ) is defined by 

^ = > (2-H) 

where the affine connection T is defined with the help of metric tensor B a p as 

rjg 7 = 2^ aS ^ 5/3 ' 7 + Bsy,f3 — Bpy,s) ■ (2.12) 

The equations Q and (|2~T0| ) can be solved iteratively, starting from a stationary point. 
In principle, the procedure to find a collective path is to find successive points at which an 
eigenvector f \ of the covariant RPA equation (|2.10| ) satisfies the force condition (|2.9| ) at the 
same time. 
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Once we get a collective path E in mult i- dimensional configuration space, a collective 
Hamiltonian is denned by evaluating the Hamiltonian (|2.5|) on the path E, 



s ^l^pl + Viq 1 ) , (2.13) 



where we assume a single collective coordinate K — 1. 

The quality of decoupling can be measured by comparing two different collective mass 
parameters that can be calculated in the theory. If we calculate the derivatives d^ a /dq l in 
terms of the tangents of the path, 

Bn - w - Bap df ■ (2 - 14) 

The other mass parameter can be obtained by using the eigenvectors f l a of the covariant 
RPA equation. 

B n = f^f^ . (2.15) 

This is equal to (-Bn) -1 if the decoupling is exact. Therefore, we define the decoupling 
measure D as 

D = [B u y l B 11 - I . (2.16) 

The size of this measure D indicates the badness of decoupling. 

If the decoupling is good, the motion orthogonal to E (motion in directions of q a ) becomes 
irrelevant in classical systems. However, it is not necessarily the case in quantum systems, 
because, according to the uncertain relation principle, we cannot require q a = (the motion 
is confined on E) and 7r a = at the same time. Therefore, one may need to include the energy 
correction into the potential V(q), which arises from quantum fluctuation with respect to 
the intrinsic degrees of freedom 0. Subsequently, the collective Hamiltonian H co i will be 
quantised in a flat space (B 11 = 1) to obtain physical quantities, such as energies and wave 
functions. 



B. Constrained local harmonic equations (CLHE) 

The most practical and straightforward way to investigate a pairing Hamiltonian would 
be to utilise the mean-field (BCS) approximation (Sec. [HI A|) in which the general prod- 
uct wave functions are no longer eigenstates of the particle number. When an intrinsic 
ground state breaks a (continuous) symmetry of Hamiltonian, generally speaking, a Nambu- 
Goldstone (NG) mode will appear. This mode corresponds to an "excitation" to connect 
degenerate vacua. In the case of particle-number breaking, this is often called "pairing 
rotation". We are interested in finding a collective coordinate orthogonal to this trivial 
mode. 

It is also well-known that, at equilibrium points, the NG mode has zero excitation energy 
and is decoupled from the other modes in the RPA order. Thus, one may easily separate 
this mode from the other physical modes (unless the physical mode happens to be zero 
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energy). However, the NG mode does not necessarily have zero energy at non-equilibrium 
points. Since we have to solve the LHE at each point on the collective path, it is not trivial 
to distinguish the NG mode from the other modes. Therefore, in this section, we provide a 
modified formulation of LHE which fixes a value of the coordinate corresponding to the NG 
mode. In the BCS parametrisation, we need the constraints to fix the particle number and 
gauge angle (Sec. HI A|) . We shall make use of the fact that the coordinate and momenta 



corresponding to the NG mode are explicitly known. 

We shall only discuss the case of a single collective coordinate K = 1 and a single NG 
mode. First we divide the set {g^} into three subsets, q 1 , q NG , and q a , a — 3, • • • ,n, which 
are assumed to be obtained by point transformations 

= f (0 , C = 9 a (q), (2.17) 

Pn = g^ a , K a = f%p fl . (2.18) 

g NG represents a coordinate for a NG mode and q 1 for a collective coordinate now we are 
trying to determine. In the case that the translational symmetry is broken, g NG and p NG 
would be the centre-of-mass coordinate and the total momentum, respectively. In the case 
of particle-number breaking, they correspond to the particle number and the gauge angle. 
In all these cases we know the explicit form of the coordinates and momenta for the NG 
mode(s). Thus, we can write g NG and p NG as functions of original coordinates £ Q and momenta 
iT a . Although in general q NG and p NG have arbitrary dependence on £ and ir, we expand them 
with respect to it up to the zero-th order for g NG and up to first order for p NG . In keeping 
with the adiabatic character of the theory, then, we get / NG (£ a ) and # ng(/ m (0)- An example 
will be shown in the next section for the particle number and the gauge angle. 

The conditions for decoupling are again given by three equations ( |2.6| ), ( |2.7| ) and 



where % = 1 and NG. In this section we use symbols ■ ■ •) representing 1 and NG. In Ref. 
|]J it is shown that the third decoupling condition ( |2.8j ) implies that the decoupled surface 
is a geodesic. The geodesic surface is defined by differential equations 

^ + r^fr3-r£-# = o, (2.19) 
fU ~ rjj/Jr + r$*/£/S = , (2.20) 

where the affine connection T is defined by Eq. fl2.12p and T in new coordinates is in the 
same manner 

f \ v - 2^ XK {P*w + ^Kv,ii - B^k) • (2-21) 

Now we consider constraints 

q NG = f G (0 = q N o G (= const.) , (2.22) 
Png = g% G K a = , (2.23) 

which freeze the NG degree of freedom. The Poisson bracket between these two constraints 
is {q ,NG ,PNc}pB = f N a9% G = I- In order to facilitate the calculation of differentiation, we 
need to define a Dirac bracket (?J by 

{F, G} DB = {F, G} PB + {F, g NG } PB {PNG, G} PB - {F, PngIpb^, G} pb . (2.24) 
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First we consider the force condition corresponding to Eq. ( |2.9| ). We differentiate the poten- 
tial V(£) = V(q) with respect to £ a keeping the constraints g NG = q^ ' and p NG = 0. 



dV 



Vif\ + VJ a a = \ ';./';, 



(2.25) 



q NG =q N G 



where the decoupling condition (|2.7| ) has been used. This can be calculated by using the 
Dirac bracket as 



dV 



<? NG =C ; 



Equating the above two equations, we obtain the force condition 

v ,a v ,ij a ■ 



(2.26) 



(2.27) 



In the case of pairing rotation, / NG may be taken as the average particle number M = (N). 
Then, V m = dV/dN can be regarded as a chemical potential fi and Eq. Q2.27p may be 
rewritten in the form similar to the constrained Hartree-Fock-Bogoliubov equation 

V a = /iAf >a + A/i . (2.28) 

In order to obtain the covariant RPA equation corresponding to Eq. Q2.10 ), we start with 
the contravariant derivative of the potential 



(2.29) 



where V' a = B^Vp and V' 1 = B 1 ^^ = B ll Vi. Again we differentiate this equation with 
respect to £ a in keeping with q NG = p NG = 0. The derivative of the l.h.s. is easily calculated 

as 



dV> 



[2.30) 



=o 



The derivative of the r.h.s. is calculated in terms of the Dirac bracket. 

<j NG =0 



/3/DB , 



(2.31) 

(2.32) 
(2.33) 



where, from the second to the third lines, we used the geodesic equation ( |2.2(J| ) and the 
covariant derivative defined by 



y,/3 — y,/3 _j_ p/3 y,1 



(2.34) 



From Eqs. Q2.30|) and ( |2.33|) , we obtain 
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{V.;! - f%giy*) /i = Vtf /i . (2.35) 

Replacing V^ 1 by lo 2 , this becomes the covariant RPA equation to determine eigenfrequencies 
uj and eigenvecotors f l a . It is worth noting that any solution f l a of Eq.( |2.35| ) with non- 
zero eigenfrequency is orthogonal to the NG mode. This may be understood as follows: 
One multiplies Eq.( |2.35| ) by and takes sum over a. Then, since the l.h.s. becomes 
V'kq — V£ G = 0, it results in 

" 2 g% G f} a = • (2.36) 

Therefore, with uj ^ 0, we have <7ng/q = 0- 

Since the NG mode emerges only when we break a symmetry of Hamiltonian, one may 
expect that the CLHE formalism should become the usual LHE when the underlying sym- 
metry is unbroken. In the next section, we will discuss the BCS approximation to derive 
a classical Hamiltonian system. In the parametrisation of Eq. ( 3.17 ) and ( 3.18Q , a state 



with = 1 and tx = corresponds to a state in a normal phase which is an eigenstate 
of particle number (from Eq.( [3.23| ), one can see the pairing gap vanish in this case). Then, 
from Eq. (|3.25|) , we have = 0. Therefore, the CLHE ( 2.27Q and (|2.35|) become equivalent 
to the LHE Q and ( pTD|) . 



The equations ( 2.27Q and ( p. 35 ) are the defining equations of CLHE formalism. As in 



the usual LHE, we require the self-consistency between these two equations: A solution of 
Eq. fl2.35|) should satisfies Eq. (|2.27| ) at the same time. In this paper, where we concentrate 



on a simple model, it is easier to solve the constraints g NG = qQ G and p NG = explicitly 
and describe everything by independent canonical variables (£*") ^a) ( a — !)■ 
However, in the realistic cases for which it is difficult to solve the constraints, this CLHE 
formalism would become more useful. 



III. TRANSCRIPTION INTO CLASSICAL FORM 

We consider a simple model where particles coupled to a harmonic oscillator interact 
through the pairing force, first introduced by Fukui, Matsuo and Matsuyanagi || in order 
to study shape mixing. This model can be regarded as a vibrating core plus valence particles 
which can move between different levels by the pairing force, 

H = H COIC + H va i , (3.1) 

where 

#core = \ (p 2 + q 2 ) , (3.2) 

#vai = E ^{q) (4c« + 4c«) - Gp]p > ( 3 - 3 ) 

a=l 

and we assume G > 0. Here a denotes the time-reversed state of a and = J2a>o c a c tt is 
the pair creation operator. The total number of levels is 2Q. The single-particle energies 
e a (q) depend on the "deformation" q of the core, which induces particle-core coupling. 



S 



We need to construct a classical Hamiltonian which describes the dynamics of the valence 
particles. In order to do this, time-dependent mean-field theory, which is known to be a form 
of Hamilton's equation, is utilised in the next section. This is of practical interest because 
it has applicability to realistic systems. For the current model one can also adopt a direct 
parametrisation of the exact wave function to construct a classical Hamiltonian (Sec. [Tl FBI) . 

A. The time-dependent mean-field (TDMF) equation 

In this section we describe a canonical parametrisation of a general product wave func- 
tion, discussed in Ref. ||. We start with a BCS state \<p ) which is the vacuum to quasipar- 
ticle operators a a , 



CL a U a C^ a V a C a , 

ot u a c^ -\- v a c a . 



(3.4) 
(3.5) 



According to Thouless' theorem [[nj, any other general product wave function \z) which is 
not orthogonal to |0o) may be written in the form 



ex p < l ^Y, z ^ a W/3 ( i <■'<■</ ■ 

Z a/3 



(3.6) 



where z a/3 = —zp a . 

The time-dependent BCS equation can be obtained by taking variations with respect to 
z and z* of the classical action, 



S 



dt 



(z 


id t - 


-H 




(z 


z) 



(3.7) 



We can introduce canonical coordinates 7r a/ g) by /3 a p = (^, a f3 + i^ a f3)/ 'v2, for a > (3. For 
a < (3 antisymmetry tells us that the same coordinates are involved, 



we can easily relate (3 and z by 



(3.8) 



0, 



[ a3 



z[l + z l z 



a/3 



(3.9) 



This means that quasi-density matrices and gnasz-pairing tensors can be expressed by 



Pa"/ 



K 



cry 











(z\z) 










(z\z) 




/3(l-/? f /?) 1/2 
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(3.10) 
(3.11) 



Using Eqs.(p.4|) and ( |3.5| ), one can relate these quantities with real density matrices p ai = 
(c^Ca) and pairing tensors n ai = (c 7 c a ). 

Then, the equations of motion take the canonical form (using an implicit matrix notation) 
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7T 



dir 

on 



with a classical Hamiltonian 



(z\H\z) 



(3.13) 



(3.14) 



Since monopole pairing interactions in nuclei were originally introduced in order to give 
a simple model of the short-range attractive force, often only the Hartree terms are taken 
into account in the mean-field calculations. Thus, we shall neglect the exchange terms in 
Sec. [IV A|. If we neglect the exchange terms for the Hamiltonian Q3.3|), only (017) = (ad) 



components are available and the matrix relations can be regarded as ordinary c-number 
relations. For instance, Eqs. (|3.10| ) and (|3 . 1 1|) now take the form 



Pc 



3 _fi - 3* 3 -W 2 



(3.15) 



where there is no summation with respect to the indices (a, a). 

The coordinates defined above depend on the reference state |0o)- F° r practical applica- 
tions, it is probably convenient to use a local coordinate at each point on a collective path 
when one solves the LHE 0. Namely, when moving to a nearby point, one changes the 
reference state to the slater determinant at this new point, making the new point the origin 
of coordinate system, and solves the LHE. However, in this paper, in order to visualise a col- 
lective path in the configuration space, we use a fixed reference state and a global coordinate 
system. Since this global coordinate system cannot describe a state \<p) which is orthogonal 
to the reference state |0o), we need to select a reference state suitable for describing all 
relevant states. 

In case that we consider a system with N = Q particles, we may take the BCS state 
with u a = v a = l/y/2 for all a as the |0o) ((^ol^l^o) = ^)- The classical Hamiltonian for 
valence particles ( |3.3| ) may be written as 



a>0 



l + ^jl-^^ + O} 



1/2- 



G 
"4 



n - e - 7T- 



1/2' 



(3.16) 



where we used notations £ a = ^ aS , tt c 



Tin 



tt 2 = Ec^, and neglected the 



exchange terms. Expanding this Hamiltonian up to second order in momentum, we obtain 
an adiabatic Hamiltonian in the form of Eq. (|2.1|) with 



d 2 H veL \ 



dlT a d7Tp 



7T = 



€q(g) €a 
a/2 ^ a 



+ G{Q 



V(£) = Wvai(e, * = 0) = y/2 J2 eMZa 



a>0 




(3.17) 
(3.18) 
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where 

'—'a 

The average number of particles is 

(z\N\z) 



1 - -£ 



1 \ 1/2 
1 -2 



Af(£,7r) 



The pairing gap is 
A(£,tt) = G 



1/2 



a>0 





p 




(* 





G 


ft 




~ ~2 




G 


ft 




2 





1/2' 



a>0 



7T 



Thus, the gauge angle in the adiabatic limit can be expressed by 



(f = arctan 



ImA 



V2 



ReA Q-ta 



7T a S a + 0(tt 3 ) 



(3.19) 

(3.20) 
(3.21) 

(3.22) 
(3.23) 

(3.24) 



>o 



Since we consider a system with N = ft particles, we may regard (A/ - — ft)/2 as g NG and 
if as p NG in Sec. |11 B| . Using the adiabatic expression of the particle number ( [3.2 lj) and the 
gauge angle (|3.24j) , we can identify 



r G (0 

dg a 



dq h 



(0 



2 


-ft 










drc a 


7T=0 





One can see that the canonical relation 



'A/ 1 



'ad 



ad 



<9/ NG 



PB ? 0& 9?* 



(3.25) 
(3.26) 

(3.27) 



is satisfied. 

In case of iV = ft = 2, the conservation of particle number in the adiabatic limit (|3.21|) 
is simply satisfied by 



6 = -6 , 

Then, 7i va i can be expressed as 

1 



7i"l = -vr 2 . 



B 



val 



42 



G{i-e 



V(0 = 2V2e£Z - G(l - £ 



2\2 



(3.28) 

(3.29) 

(3.30) 
(3.31) 
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where we assume e = e± = —62 and 

/ 1 \ 1/2 

£ = 6 = "6 , VT EE 7T! - 7T 2 = 2^ , 3 EE (l - -£ 2 J . (3.32) 

The potential is stationary at 

£ = ±1, (3-33) 

a,ld -7l(/^l ± V / ^l) • for|f|<G ' (3 ' 34) 

At £ = 1 (—1) two particles occupy the a = 1 (a = 2) level. For e ^ 0, one of the 
states corresponds to a local minimum and the other to a local maximum. For |e| < G, the 
potential has another local minimum which corresponds to a superfluid state (the upper and 
lower signs in Eq. (|3.34|) give the same state in the limit of it = 0). 



B. The time-dependent Schrodinger equation (TDSE) 

In case of iV = Q = 2, one can easily find a canonical parametrisation for the exact 
time-dependent Schrodinger equation (TDSE). With a notation \aa) = c^cjj |0), normalised 
states can be written in general as 

\iP(t)) = e im cos(a(t))|H) + e~ im sin(a(t))|22) . (3.35) 

An action, similar to Eq.( [3.7|) , may be defined by 

S = [ f dt(i)(t)\id t - H\ijj{t)) . (3.36) 
Ju 



The time-dependent Schrodinger equation can be given in terms of the least action principle. 
For the Hamiltonian if va i, we can introduce canonical coordinates (£, p)|] as follows. 

£ = a , p = —20 sin 2a . (3.37) 



The canonical equations of motion ( 3.12 ) and ( 3.12 ) hold for the Hamiltonian 



Kai = M^vailV') = 2ecos2£ - G 1 1 + sin 2£ cos f-^f) } > ( 3 - 38 ) 



where e = e± = —t2- The adiabatic Hamiltonian (|3.29|) is obtained by expanding up to 



second order in momentum. The resultant mass parameter and potential are 

BUO = , (3.39) 

V(C) =2e cos 2£-G(l + sin 2£) . (3.40) 



1 In order to avoid confusion with it = 3.14 in this section we use the symbol p for the 
momentum conjugate to £. 
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1 



q 

FIG. 1. A schematic representation of the single-particle energies as a function of the coordi- 
nate q, for the case eo = 0, x > 0- 

In this parametrisation, the mass does not depend on the single-particle energy e. It is 
apparent from Eq.( ^.35| ) that £ = (vr/2) corresponds to a pure configuration |11) ( 1 22) ) . 
There are two stationary points of the potential in a region < £ < tc which are specified 
by a condition 

G 

tan2£ = tan2a = . (3.41) 

2e 



These give eigenstates ( |3.35| ) of H va \ with = 0. For G > 0, the state in a region < £ < n/2 



always has lower energy and the mass parameter B va \ is non-negative in this region. 



IV. RESULTS 

In this section, we present the numerical results for a simple case with N = Q = 2. For 
single-particle energies, a linear dependence on the coordinate q is assumed, 

e a (q) = (-) a+1 e(q) = (-r +1 (xq-e ) , for a = 1,2. (4.1) 

The two levels cross each other at q = €q/x- A schematic figure of single-particle levels are 
displayed as functions of q in Fig. [l|. Each level has two-fold degeneracy. 

The model can be solved quantum mechanically using a direct product basis |A^ osc ) core ® 
|aa)vai where |iV osc ) core is an eigenstate of the core Hamiltonian with the oscillator quanta 
iVosc and |aa;) V ai = <4,c4|0) (a = 1,2). We use core states with iV osc < 20 for diagonalising 
the total Hamiltonian, which can be written in a matrix form 



H 



where I is the 2x2 unit matrix. Results obtained with this method will be referred below 
as "exact calculation". 



A. TDMF parametrisation 

Now we apply the theory of ALACM to this model with the mean-field parametrisation 
in Sec. [Ill A[ The Hamiltonian has the form 
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q q 

FIG. 2. Contour plots of the energy surface for the potential V(q, £), Eq.( f4.4| ). The parameters 
eo = and x = 1 are used. The left figure (a) is the energy surface for G = 0.2 and the right 
(b) for G = 1.9. Contour lines are displayed for AV = 0.5. The thick lines indicate the collective 
paths obtained by the ALACM theory. 



H = -(p 2 + B val 7r 2 )+V(q,0 , (4.3) 
V(q, = \q 2 + 2V2( X q - e )£H - G{\ - £ 2 ) 2 , (4.4) 

where B va \ is given by Eq. ( |3.3(J[ ) with e = xi ~ e o- 

Taking e = and x = 1, the potential V(q,£,) has two local minima at (q, £) = (—2, 1) 
and (2,-1) which correspond to valence-particle configurations |11) and |22), respectively. 
A collective path is now obtained by using the LHE algorithm described in Sec. O. We 
start from a local minimum, solve the RPA equations and search for the next point in the 
direction of the lowest RPA mode. Since this system has only two degrees of freedom, we 
can easily visualise the path on the 2-dimensional surface. In Fig. |2|, the potential energy 
surface and the obtained collective path are shown for two strengths of the pairing force 
G = 0.2 and G = 1.9. In the case of eo = 0, the potential landscape has a symmetry with 
respect to a rotation of 180° about the origin, V(q, £) = V(—q,—£). For a weak pairing 
force (Fig. @ (a)), each local minimum has an independent collective path which represents 
a harmonic oscillator with a fixed valence configuration (both particles in level 1 or level 2). 
These represent diabatic solutions which do not mix in this approximation. On the other 
hand, for a strong pairing force (Fig. || (b)), we get a collective path which changes the 
particle configuration and connects the two local minima. 

We have thus found that the system automatically selects different type of collective 
paths according to the strength of pairing force. In this approximation, we get completely 
diabatic solutions for G = 0.2 which will result in no parity splitting upon quantisation. Of 
course, for G ^ 0, an exact calculation always gives some splitting in energy. However, the 
splitting turns out to be extremely small for the weak pairing case (see below). 

Next, in order to discuss a case of asymmetric potential landscape, let us take e = 0.4, 
X = 0.7 and G = 0.1. Fig. |] (a) shows the potential energy surface and collective path. The 
potential, which is no longer symmetric about the 180° rotation about the origin, has a deep 
minimum at ((?,£) = (—1-4, 1) and a shallow minimum at q — (1.4, —1). Each minimum has 
again an independent collective path parallel to the horizontal line (g-axis), representing the 
diabatic dynamics. In Fig. |3] (b), we show results of exact calculation for this parameter 
set, together with the diabatic potential curves corresponding to the valence configuration 
1 11) and |22). The number on each level indicates the percentage of the main configuration 
in the wave function. The exact calculation shows strong diabaticity of the system and two 
different configurations hardly mix. This agrees with the independent paths obtained in the 
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FIG. 3. (a) The same as Fig. | but for e = 0.4, x = 0.7 and G = 0.1. (b) Diabatic potential 
curves as functions of q and eigenenergies of the exact calculation. The number on each level 
indicates the percentage of the main configuration in the wave function. 




FIG. 4. Born-Oppenheimer potentials V(q) (solid lines), Eq. (|4.6| ), for G = 0.2 (a) and for 
G = 2 (b). The parameters eo = and x = 1 are used. The dotted lines indicate the potential 
including the scalar gauge potential, V(q) + &(q). See text for details. 

ALACM. Therefore, one may state that the ALACM combined with the BCS approximation 
can reproduce the qualitative feature of the pairing systems. 



B. Born-Oppenheimer (BO) approximation 

Before discussing another parametrisation in the next section, let us review the Born- 
Oppenheimer (BO) approximation. The formulation in the next section can be seen as the 
generalised version of BO theory, in which a collective coordinate is chosen self-consistently. 

Using the conventional adiabatic theory which assumes motion of the valence particles 
is much faster than motion of the core, the effective Hamiltonian can be obtained by diago- 
nalising H va _\ at each value of q. 

HBO = ^ + V BO {q) , (4.5) 

2 

Vbo{q) = y - V^xg-eo^ + G 2 - G . (4.6) 
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If we neglect the interaction (G = 0), the potential Vbo(q) has two local minima at q — ±x 
corresponding to different valence-particle configurations. Fig. |] shows the BO potentials 
Vbo(?) f° r weak (Fig. | (a)) and strong pairing force (Fig. [| (b)) with x = 1 and e = 0. 
In Ref. |H it has been shown that this approximation works well for the strong pairing case 
but it considerably overestimates the mixing between two minima for the weak pairing case. 

One possible way to improve the BO theory is to determine a collective path self- 
consistently and take into account the effect of mass parameter B va \. This will be discussed 
in the next section. The other, which is discussed here, is to take account of the (Berry's) 
"gauge" potentials which emerge from the derivative of the eigenstates of the fast degrees 
of freedom (valence particles) \n) with respect to the slow coordinate q. The new effective 
Hamiltonian will be 

Hbo = ~ (P - A n (q)f + V BO (q) + $ n {q) , (4.7) 

where 

A n {q) = t(n\d g n) , (4.8) 
$ n (q) = l(d q n\(l-\n)(n\)\d q n) . (4.9) 

In the present case, we take |n) to be the lowest eigenstate of the valence Hamiltonian which 
can be expressed as in Eq. (|3.35| ), with a given by Eq.( j3.41f) . Then, the gauge potentials are 



calculated as 

A(q) = 0, iXG) * rn2 - (4-10) 

2{4(xg-eo) 2 + G 2 } 

This shows that there is no Berry's phase for the current problem, but does imply an addi- 
tional scalar potential. The additional potential is positive definite and has the maximum 
value at the crossing point q = e /x- This scalar potential is shown in Fig. f| by dotted 
lines. The peak height at a level crossing is getting higher for weaker pairing. This term 
will suppress the mixing between two minima for the case of weak pairing, and improve the 
result of quantisation (see below). 



C. TDSE parametrisation 



In the Sec. |IV A| , we have shown the results using the mean-field parametrisation of 
the valence Hamiltonian, which leads to a diabatic solution in the weak pairing case. In 
this section, we use the parametrisation described in Sec. |II [B[ This conserves number of 
particles exactly and can be regarded as a generalised BO theory. 

The total Hamiltonian is in a form 

H = i(p 2 + iW)+V(g,0, (4- 11 ) 
V(q,0 = lq 2 + 2( X q-e )cos2^-G(l + sm20 , (4.12) 

where S va j is given by Eq.( |3.39| ). 
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FIG. 5. Contour plots of the potential energy V(q,£), Eq.(p~T2|), for G = 0.2 (a) and G = 2 
(b). The parameters eo = and x = 1 are used. The thick lines indicate the collective paths 
obtained by the ALACM, while the dotted lines are the paths in the BO approximation. 



With €q = and x = 1> the potential surface has again two local minima at 

' y/16 - G 2 1 fVlQ-G r 



, - arccos 



4 



, [VIQ-G 2 1 / VIQ-G r 
and I , - arccos I 



(4.13) 



and a saddle point at (0, 7r/4) where we restrict £ to the region < £ < it/2. In this 
parametrisation, £ represents the mixing angle between the two valence configurations, |ll) 
and 1 22). If we minimise the potential with respect to the coordinate £ at each value of q, 
the potential V(q, £o(?)) becomes the same as the BO potential (|4.6|) . If G = 0, the local 
minima correspond to the pure configurations, |11) (£ = 0) and 1 22) (£ = tt/2). In general, 
these minima have finite values of mixing angle, which was not the case in the mean-field 
parametrisation. 

Fig. |5| shows the potential energy surfaces and the 1-dimensional collective paths. For a 
strong pairing force (Fig. [| (b)), we have a smooth path which connects two local minima. 
On the other hand, for the case of weak pairing (Fig. |5] (a)), the path exhibits a peculiar kink- 
like behaviour near q = ±0.8. Close to a local minimum point, the collective path represents 
a harmonic oscillator (the core Hamiltonian) which is almost parallel to a g-axis (horizontal 
line). On the other hand, at the saddle point (<?,£) = (0,7r/4), the collective path is trying 
to go through the valley of potential. These paths are not able to be connected smoothly 
for the weak pairing case, which results in the strange back-bending. In the vicinity of this 
back-bending, the decoupling from the other degree of freedom is very bad (see discussion 
below) . 

In the same figure, we show a path (dotted line) which is obtained by minimising the 
potential energy at each value of q (BO approximation). For weak pairing, the path obtained 
by the ALACM theory is found to be quite different from that by the BO theory. On the 
other hand, for strong pairing, we can hardly see the difference between them. In the BO 
approximation, the motion along £-axis is assumed to be much faster than that along g-axis. 
In other words, the mass of valence Hamiltonian 1/B va \ should be much smaller than the 
mass of core Hamiltonian, which can be expressed as the condition 



G > sin2£ . 



(4.14) 



It is obvious from this equation that the BO approximation fails for the weak pairing force 



17 




FIG. 6. Potentials along the paths (thick lines) displayed in Fig. |5| as functions of a collective 
coordinate s for G = 0.2 (b) and G = 2 (d). Solid (dotted) lines are the potentials without (with) a 
scalar gauge potential <3?(s). The upper part, (a) and (c), shows the decoupling measure D defined 
in Eq.( 2.16p for G = 0.2 and G = 2, respectively. 



In order to requantise the obtained collective Hamiltonian, it is convenient to define the 
collective coordinate s on the path so as to make the collective mass parameter B 11 = 1. 
We use the relation 

to renormalise the effect of collective mass into a new coordinate s. The requantised collective 
Hamiltonian will be in a form 

n = + (4.i6) 

where V(s) = V(q(s),£(s)). 

In Fig. [], the potential and the decoupling measure D ( f2.16| ) are shown as functions of 
coordinate s. Comparing these potentials with the BO potentials in Fig. |[ for the strong 
pairing case, we see that they resemble each other and the decoupling is found to be very 
good everywhere on the path. For the weak pairing force, however, the distance between 
two minima of the potential is larger in the ALACM than in the BO theory. The decoupling 
is extremely bad where the path has a kink (the decoupling measure D is order of 10 3 at 
the peak), though the decoupling is relatively good anywhere else on the path. 

As we have done for the BO approximation, we can incorporate the effective "gauge" 
potentials into the ALACM theory. We assume gaussian wave functions for the intrinsic 
(non- collective) degrees of freedom ("harmonic approximation"). Here, collective and intrin- 
sic coordinates are expressed by s and x, respectively. The wave function is approximated 
in a form of semi-direct product, 

*(s,ar) = Vccrt(s) ® ^mtr(s, x) , (4.17) 
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where the wave functions for the intrinsic motion if)int T (s,x) are assumed to be a gaussian, 

Ms , x)= (^y\ xv ^^M x ^ . ,4.18) 

Here the frequency u x (s) are calculated by the covariant RPA at each point on the collective 
path. Since the lowest RPA mode is supposed to be along the collective path, u x (s) is the 
frequency of the second RPA mode. 

Identifying ^mtrts, x) with a state \n) in Eqs.( |4.8|) and (|4.9|), the (Berry's) gauge poten- 
tials become 



<9V>infa 



A(s) = i dx^ t ^L = , (4.19) 



1 r°° 
$(s) = - / (hi 



intr ds 

2 

^intr 



oo 



ds 



±(ll„. W ) . (4.20) 



Again, the vector potential A vanishes and the scalar potential is positive definite. In Fig. |6|, 
we show this scalar gauge potential as a function of s (dotted lines). For a strong pairing 
force $(s) ~ everywhere (one cannot see the dotted line in the figure because it almost 
coincides with the solid line), however, for weak pairing, it has sharp peaks at kinks where 
the decoupling measure D also has large values. The requantised Hamiltonian (|4.16|) will 
be modified into 

K = -l-^ + V(s)+m- (4.21) 



The wave function ip co \(s) in Eq. (|4.17|) is required to be an eigenfunction of this modified 
Hamiltonian. 

The results after requantisation are summarised in Fig. [7[ There we show the spectra of 

(1) the BO theory without gauge potentials, (2) the BO theory with gauge potentials, (3) 
the ALACM theory without gauge potentials, (4) the ALACM theory with gauge potentials, 
and (5) the exact calculation. For a strong pairing force, all approximation schemes work 
well. On the other hand, for a weak pairing force, the result (1) overestimates the splitting 
of parity doublets. This overestimate is corrected by including the scalar gauge potential 

(2) . The same effect also appears for ALACM results (3) and (4). However, in this case, the 
suppression of parity splitting originates in the mass parameter, not in the gauge potentials. 



Since we renormalise the mass parameter with Eq. (|4.15|) , this produces a larger distance 
between two local minima of the potential V(s), as compared to Vbo(<?)- The tunnelling 
probability between two minima is suppressed by this mass-renormalisation effect. The 
effect of gauge potentials, the difference between (3) and (4), turns out not to be significant 
in the ALACM. The agreement with the exact calculation is good for low-lying states (up 
to the third excited state), though, as expected, it is not so good for higher- lying states. 
We also show transition matrix elements | (n'\q\n) | in the same figure. They are calculated 

by 

/oo 
dq ip*>(q)qM<l) for BO, (4.22) 
-oo 

/■OO 

(n'\q\n) = ds </(s)g(s)^ n (s) for ALACM , (4.23) 
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FIG. 7. Spectra up to the fifth excited states for (1) the BO theory without &(q), (2) the BO 
theory with $(g), (3) the ALACM theory without $(s), (4) the ALACM theory with 3>(s), and 
(5) the exact calculation. Figure (a) is for G = 0.2 and (b) for G = 2. The numbers next to the 
arrows denote The values of the transition matrix element |(ra|g|n + 1)|. 

where the wave functions ip(q) and ip(s) are normalised with respect to the coordinates q 
and s, respectively. Again, for the case of weak pairing, the amplitude between the first 
and second exited states is overestimated in BO theory (1) due to too strong mixing of two 
minima. This is corrected in ALACM (3), (4) and in BO with gauge potential (2). 

Table | shows the expectation values of the decoupling measure D with respect to indi- 
vidual states, 

/oo 
ds iP*(s)D(s)iP n (s) . (4.24) 
-oo 

These values are about 7 orders of magnitude smaller for the case of strong pairing G = 2 
compared to the weak pairing G = 0.2. For the weak pairing force, it turns out that all 
states considerably suffer bad decoupling. Surprisingly, still we have obtained reasonable 
agreement for the low-lying spectra and transition amplitudes. 

V. CONCLUSION 

We have investigated a simple system with level crossings where the monopole pair- 
ing interaction plays an important role and the potential landscape exhibits multiple local 
minima. To apply the ALACM theory, canonical coordinates and a classical Hamilton's 
equation have been introduced by parametrisation of the time- dependent states. We have 
done this in two different ways. One is based on the time-dependent mean-field (BCS) theory 
in which states are not eigenstates of particle number. It is inevitable that this introduces 
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TABLE I. Expectation values of the decoupling measure (n\D\n) with respect to individual 
eigenfunctions for cases (3) and (4) in Fig. [?]. See the main text and the caption of Fig. [?] for more 
details. 



n 


(3) 


G = 0.2 


(4) 


(3) 


G = 2 


(4) 


1 


57 




45 


3.97 x 10" 6 




3.98 x 10~ 6 


2 


55 




45 


1.72 x 1(T 6 




1.71 x 10~ 6 


3 


172 




162 


3.17 x 1(T 6 




3.17 x 10~ 6 


4 


181 




167 


3.19 x 1(T 6 




3.19 x 10~ 6 


5 


8.0 




7.1 


1.01 x 10~ 6 




1.00 x 10~ 6 


6 


119 




125 


2.99 x 10~ 6 




2.99 x 10~ 6 



a Nambu-Goldstone (NG) mode (pairing rotation) in the theory. We have developed the 
constrained local harmonic formalism which can treat this spurious component and find a 
path orthogonal to the NG mode in the configuration space. The other is based on the 
exact time-dependent Schrodinger equation. It is possible to derive a classical Hamilton's 
equation for a simple two-level system by regarding the mixing angle and the relative phase 
as canonical variables, though this becomes more complicated when the number of levels 
increases. In this parametrisation, states are always eigenstates of particle number and there 
is no NG mode. 

Using the mean- field parametrisation (neglecting Fock terms), we have found that the 
system exhibits completely different collective paths depending on the strength of pairing 
force. For a weak pairing force, the collective path cannot connect one local minimum with 
another. On the other hand, a path which connects two local minima is obtained for the 
case of a strong pairing force. In this model, different local minima correspond to different 
configurations of valence particles. Therefore, the strength of pairing force determines either 
adiabatic or diabatic behaviours of the system. 

Using the parametrisation of the exact time-dependent Schrodinger equation, if G ^ 0, 
a collective path starting from a local minimum always passes through a saddle point and 
reach the other minimum. For the strong paring force, the path is very smooth and is well 
decoupled from the other degree of freedom orthogonal to the path. However, for the weak 
pairing force, the path shows a peculiar behaviour (back-bending) on the way where the 
decoupling is very bad. This is qualitatively consistent with the mean-field results. 

We have requantised the 1-dimensional collective Hamiltonian and compared the spectra 
and transition matrix elements with exact calculations. For the strong pairing force, the 
agreement is excellent for any kind of adiabatic approximation (BO, ALACM). However, for 
the weak pairing, the BO approximation fails to reproduce the qualitative features of exact 
spectra, while the ALACM well reproduces the low-lying spectra and transition amplitudes. 
The agreement becomes worse for the higher-lying states beyond a barrier height between 
two local minima. The theory is based on the adiabatic assumption in which Hamiltonians 
are expanded with respect to momentum up to second order. Thus, this disagreement may 
come from the higher-order terms we have neglected, because the higher-lying states are 
supposed to have larger momenta. 

We have also discussed the effect of gauge potentials arising from change of intrinsic 
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states along the path. Inclusion of these potentials turns out to improve the BO calculations 
considerably but to leave the ALACM results almost unchanged. 

The authors of Ref. || have studied the same model Hamiltonian and concluded that 
for the case of a weak pairing force the adiabatic approach (BO approximation) is not 
appropriate for describing the system and the diabatic approach is needed. However, one 
should bear in mind that, in the BO approximation, the collective coordinate q was chosen 
by assuming that the motion of valence particles described by (£, 7r) is much faster than the 
motion along q. This is apparently not the case for the weak paring G<C 1. In our adiabatic 
approach, all degrees of freedom are treated equivalently, taking account of the mass tensor 
in the multi-dimensional configuration space. Thus the collective path is very much different 
from the one in the BO approximation for the case of weak pairing. Actually, the collective 
path turns out to represent the completely diabatic dynamics in some cases. Therefore, if 
one has a decoupled path determined by the system itself, the adiabatic theory may account 
for the diabatic dynamics. 

For the practical applications of the ALACM theory, generally speaking, it is not easy 
to solve the constraints about the NG mode explicitly. Since nuclei show many kinds of 
symmetry breaking (translational symmetry, rotational symmetry, gauge symmetry), the 
constrained local harmonic formalism discussed in Sec. [11 B| will be very useful in such cases. 
We shall apply this method to more realistic problems in the near future. 

In conclusion, the ALACM improves on the conventional BO theory. It works well even 
when the mass of valence Hamiltonian is smaller than that of the core Hamiltonian, because, 
in principle, the ALACM treats all degrees of freedom equivalently. The ALACM fails to 
reproduce properties of higher-lying states in the case of weak pairing. This is definitely due 
to a limitation of the adiabatic theory. Upon requantisation, a reasonable agreement with 
the exact calculation has been achieved in both cases of strong and weak pairing forces, at 
least for the low-lying states. 



ACKNOWLEDGMENTS 

This work was supported by a research grant from the Engineering and Physical Sciences 
Research Council (EPSRC) of Great Britain. 



22 



REFERENCES 



[1] A. Klein, N.R. Walet and G. Do Dang, Ann. Phys. 208 (1991) 90. 
[2] N.R. Walet, G. Do Dang and A. Klein, Phys. Rev. C43 (1991) 2254. 
[3] N.R. Walet, A. Klein and G. Do Dang, Phys. Rev. C45 (1992) 249. 
[4] D.L. Hill and J.A. Wheeler, Phys. Rev. 89 (1953) 1102. 
[5] R. Bengtsson and W. Nazarewicz, Z. Phys. A334 (1989) 269. 

[6] G. Bertsch, Frontiers and borderlines in many-particle physics, (Enrico Fermi School, 
CIV Corso, 1988), p.41. 

[7] P.A.M. Dirac, Lectures on Quantum Mechanics, (Belfer Graduate School of Science, 
Yeshiva University, New York, 1964). 

[8] T. Fukui, M. Matsuo, K. Matsuyanagi, Prog. Theor. Phys. 85 (1991) 281. 

[9] J. -P. Blaizot and G. Ripka, Quantum Theory of Finite Systems, (The MIT Press, Cam- 
bridge, Massachusetts, 1986), p. 298. 
[10] D.J. Thouless, Nucl. Phys. 21 (1960) 225. 



23 



